Detailed Jarzynski Equality applied to a Logically Irreversible Procedure 
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A single bit memory system is made with a brownian particle held by an optical tweezer in a 
double-well potential and the work necessary to erase the memory is measured. We show that the 
minimum of this work is close to the Landauer's bound only for very slow erasure procedure. Instead 
a detailed Jarzynski equality allows us to retrieve the Landauer's bound independently on the speed 
of this erasure procedure. For the two separated subprocesses, i.e. the transition from state 1 to 
state and the transition from state to state 0, the Jarzynski equality does not hold but the 
generalized version links the work done on the system to the probability that it returns to its initial 
state under the time-reversed procedure. 
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The connection between thermodynamics and infor- 
mation is nowadays a widely studied problem [1-5]. The 
main questions concern the amount of energy necessary 
in order to perform a logical operation in a given time and 
how the information entropy is related to the free energy 
difference between the initial and final state of this log- 
ical operation. In this context the Landauer's principle 
[6] is very important as it states that for any irreversible 
logical operation the minimum amount of entropy pro- 
duction is — fcs ln(2) per bit commuted by the logical op- 
eration, with ks the Boltzmann constant. Specifically a 
logically irreversible operation is an operation for which 
the knowledge of the output does not allow to retrieve 
the initial state, examples are logical AND, OR and era- 
sure. In a recent paper [8] we have experimentally shown 
that indeed the minimum amount of work necessary to 
erase a bit is actually associated with this Landauer's 
bound which can be asymptotically reached for adiabatic 
transformations. The question that arises naturally is 
whether this work corresponds to the free energy differ- 
ence between the initial and final state of the system. 
To answer to this question it seems natural to use the 
Jarzinsky equality [7] which allows one to compute the 
free energy difference between two states of a system, 
in contact with a heat bath at temperature T. When 
such a system is driven from an equilibrium state A to a 
state B through any continuous procedure, the Jarzyn- 
ski equality links the stochastic work Wst received by the 
system during the procedure to the free energy difference 
AF = Fb — Fa between the two states: 

Where (.) denotes the ensemble average over all possible 
trajectories, and /3 = (see eq. 2 for the precise 

definition of the work Wst)- 

In this letter we analyze the question of the applica- 
tion of eq. 1 for estimating the AF corresponding to 
the erasure operation in our experiment, in which a col- 
loidal particle confined in a double well potential is used 
as a single bit memory. We will show that a detailed 



Jarzynski Equality is verified, retrieving the Landauer 
limit independently of the work done on the system. 

The setup has already been described in a previous 
article [8] and we recall here only the main features. 

A custom-built vertical optical tweezers is used to re- 
alize a two-state system: a silica bead (radius R = Ifim) 
is trapped at the focus of a laser beam (wavelength 
lQ24nm) which is rapidly switched (at a rate of lOkHz) 
between two positions (separated by lAbfim) using an 
acousto-optic deflector. A disk-shaped cell {18mm in di- 
ameter, 1mm in depth) is filled with a solution of beads 
dispersed in bidistilled water at low concentration. The 
bead used for the experiment is trapped by the laser and 
moved into the center of the cell (with gap ^ 80/im) 
to avoid all interactions with other beads. The trapped 
bead feels a double-well potential with a central barrier 
varying from 2kBT to more than SksT depending on the 
power of the laser (see figure 1, a and b). The left well 
is called "0" and the right well "1" . The position of the 
bead is tracked using a fast camera with a resolution of 
108nm per pixel, which after treatment gives the position 
with a precision greater than lOnm. The trajectories of 
the bead are sampled at 502Hz. 

The logical operation performed by our experiment is 
the erasure procedure. This procedure brings the system 
initially in one random state (0 or 1 with same probabil- 
ity) in one chosen state (e.g. 0). It is done experimentally 
in the following way. 

At the beginning the laser power is high {48mW) so that 
the central barrier is more than 8kT and the characteris- 
tic jumping time (Kramers Time) is about 3000s, which 
is long compared to the time of the experiment, and the 
bead is trapped in one well-defined state (0 or 1). The 
laser power is first lowered (in a time Tioweriiig — Is) 
to 15mW so that the barrier is about 2.2kBT and the 
jumping time falls to 10s. A viscous drag force linear 
in time is induced by displacing the cell with respect to 
the laser using a closed-loop 3-axis NanoMax stage. The 
force is given by / = —jv where 7 = GirRrj (rj is the vis- 
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FIG. 1: Measurement of the potential in which the sihca bead 
is confined, with no external force for two different laser pow- 
ers (a: 48mW and b: 15mW), and extrapolation for low 
power with two different value of the force (c: 1 x 10~^'''N 
and d: 3 X lO'^'^A'^). The potential is computed from the 
equilibrium probability density function of the bead's posi- 
tion. 



cosity of water) and v the speed of displacement. It tilts 
the double-well potential so that the bead ends always 
in the same well (e.g. state 0) independently of where it 
started. At the end, the force is stopped and the central 
barrier is raised again to its maximal value (in a time 
Trising = Is). Between two successive procedures the 
system is left to equilibrate for 4 s. The experimental 
procedure is sketched in figure 2. A procedure is fully 
characterized by its duration t and the maximum value 
of the force applied fmax- Its efficiency is characterized 
by the "proportion of success" Ps, which is the propor- 
tion of trajectories where the bead ends in the chosen 
well (e.g. 0), independently of where it started. 
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FIG. 2: Schematic representation of the erasure procedure. 

The ideal erasure procedure is a logically irreversible 
operation because the final state gives no information 



about the initial state. For one bit of memory [6], it 
corresponds to a change in the entropy of the system 
AS* — — fcBln(2). The procedure can arbitrarily be de- 
composed in two kinds of sub-procedures: one where the 
bead starts in one well and ends in the other (e.g. 1 — >■ 0) 
and one where the bead is initially in the same well where 
it should be at the end of the procedure (e.g. — > 0). 

The two accessible quantifies are x{t)^ which is mea- 
sured, and f{t) which is imposed by the displacement 
of the cell. The derivatives are estimated using the dis- 
cretization x{t + At/2) w [x{t + At) ~ x{t)]/At. Start- 
ing from these quantities it is possible to measure the 
stochastic work Wst done during the erasure procedure. 

For a colloidal particle confined to one spatial dimen- 
sion and submitted to a conservative potential V{x,X), 
where A — X(t) is a time-dependent external parameter, 
one can define the stochastic work received by the system 
along a single trajectory [5]: 

Wst[x{t')] = I^Adi (2) 

Here the potential is V{x,t) — Uo{x,I{t)) — x x f{t), 
where Uo is due to the optical trapping and /(<) is the 
intensity of the laser (see figure 1). If the bead does not 
jump from one well to the other during the modulation 
of the height of the barrier this part of the procedure 
does not contribute to the work received by the bead. 
Then the work can be computed only when the external 
force is applied (between t — and t — t). The force 
is directly the control parameter, and considering that 
f{t = 0) = = f{t = r), it follows that the stochastic 
work is equal to the classical work W: 

Wst[x{T)]^ f -xfdt^ f fxdt^W[x{T)] (3) 
Jo Jo 

The two integrals have been calculated for all the tra- 
jectories of all the procedures tested. Amongst all of 
them, the maximal difference \ Wst — W\ observed was of 
O-OGhsT, which is negligible. 

We now analyse the results of our experiments. For 
every chosen duration r, the maximal force fmax was set 
to different values (typically between 1 and 6 x 10"^"*^^). 
For each set of parameters {T,fmax), the procedure was 
repeated several hundreds of times to be able to com- 
pute statistical values. For each r, the value of fmax is 
optimized in order to be as small as possible and give a 
proportion of success Ps > 90%. 

The trajectories where the information is erased, i.e. 
the ones where the bead ends where it was supposed 
to be (e.g. in state 0), are selected. The mean of the 
work received {W)^q and the logarithm of the mean of 
its exponential — In ((e~^^)_^Q) are calculated, where 
(•)_j.o stands for the mean on all the trajectories end- 
ing in 0. We call the value - In ((e"'^'^)^;,) the ef- 
fective free energy difference AF^ff. The error bars 
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on this value are estimated by computing the mean on 
the data set with 10% of the points randomly excluded, 
and taking the maximal difference in the mean value 
observed by repeating this operation 1000 times. The 
results are shown in figure 3. The mean work {W)^q 
decreases when the duration of the procedure increases. 
For the optimized values of the force, it follows a law 
{W{t))^q = kBT\n{2) + B/t where S is a constant, 
which is the behavior for the theoretical optimal proce- 
dure [9]. A least mean square fit gives B = 8.45 fc^T.s. 
The effective free energy difference AFg/ / is always close 
to the Landauer limit kBTh\{2)^ independently of the 
value of the maximal force or the procedure duration. 
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FIG. 3: Mean of the work (*) and effective free energy differ- 
ence (x) for different procedures. The over-forced procedures 
(red) have a proportion of success Ps ~ 95%, the optimized 
procedures (black) have Ps > 91%, the under- forced proce- 
dures have Ps > 83% (except the last point, that has Ps « 
75%). The fit (blue Une) is: (W(t))^o = A;srin(2) + B/t 
with B a constant. 



The mean of the exponential can be computed on the 
sub-procedures by sorting trajectories in function of the 
initial position of the bead. Specifically : 
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where 



and Moo = (e-'^'^), 



(4) 



(5) 



The values Miq and Moo a-re plotted in fig. 4. The sum 
MiQ + Moo is always close to 1, which corresponds to the 
fact that AFeff is close to ksT ln{2), but Afoo decreases 
with T whereas Miq increases consequently. 

These results can be understood in the following way: 
Since the memory erasure procedure is made in a cyclic 
way and AS* = —ks ln(2) it is natural to await AF — 
ksT ln{2). But the AF that appears in the Jarsyn- 
ski equality is the difference between the free energy 
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FIG. 4: Mean of the exponential of the work, for the sub- 
procedures 1—^0 (blue) and — > (red). For readability 
questions, only one value is shown for each t: it corresponds 
to the procedure with the highest proportion of success Ps- 



of the system in the initial state (which is at equilib- 
rium) and the equilibrium state corresponding to the fi- 
nal value of the control parameter: F(A(r)) — F(A(0)). 
Since the height of the barrier is always finite there is 
no change in the equilibrium free energy of the system 
between the beginning and the end of the procedure. 
Then AF = 0, which implies (e~'^'^=*) = 1. Neverthe- 
less Vaikuntanathan and Jarzysnki [10] have shown that 
when there is a difference between the actual state of the 
system (described by the phase-space density pt), and 
the equilibrium state (described by p^"^), the Jarzynski 
equality should be modified: 



(x,t) 



Where is the mean on all the trajectories that pass 

through X at t. 

In our procedure, the selection of the trajectories where 
the information is actually erased, corresponds to fix x 
to the chosen final well (e.g. state 0) at the time t — t. 
It follows that p{0,t) is directly Ps, the proportion of 
success of the procedure, and /9^''(0, A(t)) = 1/2. Then: 



'I3W{t) 



1/2 
Ps 



(7) 



Similarly for the trajectories that end the procedure in 
the wrong well (e.g. state 1) we have: 



-l3W{r} 



1/2 
l-Ps 



Taking into account the Jensen's inequality, i.e. (e 



we find that equations 7 and 8 imply: 

W^o>fcBT[ln(2)+ln(Ps)] 

> fcsT [ln(2) -t- ln(l - Fs)] 



(8) 



> 



(9) 
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Notice that the mean work dissipated to reaHze the pro- 
cedure is simply: 



{W) =Psx (W)^^ + (1 - Ps) X (W)^ 



(10) 



where (.) is the mean on all trajectories. Then using the 
previous inequalities it follows: 

{W) > ksT [ln(2) + Ps In(Fs) + (1 - Ps) ln(l - Ps)] 

(11) 

which is indeed the generalization of the Landauer's limit 
for Ps < 1. In the ideal case where Ps — I, all the 
trajectories end in state at time r, and finally: 



( 



e'^^) = (e-^'^)^ =1/2. 



(12) 



Since this result remains approximatively verified for pro- 
portions of success close enough to 100%, it explains why 
in the experiment we find AF^jf ~ ksT ln{2). 

To understand the evolution of Miq and Mqq, we need 
to consider the subproccdurcs 1 — > and — > sepa- 
rately. In this case the Jarzynski equality does not hold 
because the initial conditions are not correctly tested 
(selecting trajectories by their initial condition intro- 
duces a bias in the initial equilibrium distribution). But 
Kawai and coworkers [11] have shown that for a parti- 
tion of the phase-space into non-overlapping subsets Xj 
(j — 1, K) there is a detailed Jarzynski Equality : 



Pi ' Pj 



(13) 



with: 



Pj 



p(ta) dxdp and pj 



Xj 



p{ta) dxdp (14) 



where p{ta) and p{ta) are the phase-space densities of 
the system measured at the same intermediate but oth- 
erwise arbitrary point in time, in the forward and back- 
ward protocol, respectively. Here, there are only two 
subsets i = {00,10}, defined by the position where the 
bead starts. By taking ta = the starting point of the 
procedure, we have poo = 1/2 = pio, and poo (resp. 
pio) identifies with the probability Pqo (resp. Pqi) that 
the system returns into its initial state, i.e. state 
(resp. state 1), under the time-reversed procedure. Since 
g-^AFe// — 1/2 it follows from eq. 13 and the definition 
in eq. 5 that: 
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(15) 



This result is similar to the one reported in ref. [12, 13] 
for procedures with feedback. It should be noticed that 
here Pqi + ^oo = 1- It is reasonable to think that for 
time-reversed procedures (that always start in state 0) 
the probability of returning to state 1 is small for fast pro- 
cedures and increases by increasing the duration r, which 



explains qualitatively the behavior of Mqo and Mio ob- 
served experimentally. To be more quantitative one has 
to measure Pqi and Pqo, but the time-reversed procedure 
cannot be realized experimentally, because it starts with 
a very fast rising of the force, which cannot be reached 
in our experiment. 

Thus, in order to verify eq. 15, we performed a nu- 
merical simulation, where it is possible to realize the cor- 
responding time-reversed procedure and to compute Pqi 
and Pqo. Our experimental system can be described by 
the over-damped Langevin equation: 



dV 
dx 



(16) 



where ^ is a gaussian white noise with zero mean and 
correlation (^(t)C(i')> = '2-fkBTS{t - t'). 

Simple numerical simulations were made by integrat- 
ing this equation with Euler method, for a set of proce- 
dures as close as possible to the experimental ones. Some 
results are showed in the following table: 



t{s) 


fmax (lO-l^iV) 




Poi 


Moo 


Poo 


success (%) 


5 


3.77 


0.19 


0.19 


0.84 


0.81 


97 


10 


2.83 


0.30 


0.30 


0.73 


0.70 


96.5 


20 


1.89 


0.45 


0.41 


0.63 


0.59 


94 


30 


1.89 


0.45 


0.44 


0.60 


0.56 


94.5 



The agreement between Mio (resp. Mqo) and Pqi 
(resp. Poo) is quantitative, and we also retrieve the fact 
that Mio -I- Mqo is always close to 1 for any set of parame- 
ters with reasonnable success rate, as in the experiments. 
It was also verified that for proportions of success 
< 100%, if one takes all the trajectories, and not only 
the ones where the bead ends in the state 0, the classical 
Jarzynski equality is verified: ^e"*^*^'*) = 1 (for these 
specific simulations, Tiowering and Trising were taken 
equal to 0.1s to avoid problems when the bead jumps dur- 
ing this phase of the procedure). This result means that 
the small fraction of trajectories (sometimes < 1%) where 
the bead ends the erasure procedure where it shouldn't 
is enough to retrieve the fact that AP = 0. 

As a conclusion, it has been experimentally shown 
that for a memory erasure procedure of a one bit sys- 
tem, which is a logically irreversible operation, a detailed 
Jarzynski equality allows us to retrieve the Landauer's 
bound for the work done on the system independently 
on the speed in which the memory erasure procedure is 
performed. Furthermore we show that the division of 
the procedure into two sub-procedures is useful in or- 
der to link the work done on the system to the proba- 
bility that the memory returns to its initial state under 
the time-reversed procedure. These results are important 
because they clarify the use of the Jarzinsky equality in 
irreversible operations. 
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